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Abstract 

A new QCD calculation of the mass of the nucleon is presented. It 
makes use of a polynomial kernel in the dispersion integrals tailored 
to practically eliminate the contribution of the unknown l/2 + and 
l/2 _ continuum. This approach avoids the arbitrariness and other 
drawbacks attached to the Borel kernel used in previous sum rules 
calculations. Our method yields stable results for the nucleon mass 
and coupling for standard values of the condensates. The prediction 
of the nucleon mass m^r = (0.945± .045)GeU is in good agreement 
with experiment. 
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1 Introduction 



The QCD sum rules method introduced by Shifman et al. [T] has extended 
the applicability of QCD far beyond simple perturbation theory. The method 
was adapted to the case of nucleons by Ioffe [2] and independently by Chung, 
Dosch, Kremer and Schall [3]. These authors showed how to approach one of 
the fundamental problems of QCD, the calculation of baryon masses from the 
Lagrangian and the vacuum condensates. 

Nucleon mass sum rules start with the correlation function 

n{q) = i J d 4 xe iqx (0|T7(ar)»y(0)|0) (1) 

where r\ is a nucleon interpolating held constructed from local QCD operators 
with the quantum numbers of the nucleon. We will choose [2] 

abc i a/~i A b\ A ic 

j] = e (u 07 u )7s7 a . 

which couple maximally to the nucleon. The correlator can be decomposed in 
terms of invariants, 

n(<?) = ^7 Ai n 1 (< z 2 ) + n 2 ( g 2 ) 

with 7^ standing for the Dirac matrices. TL(t — q 2 ) is an analytic function in 
the complex i-plane with a pole at t = m 2 N and a cut along the positive real 
axis starting at t = (toat + m w ) 2 . The sum rule methods can be traced back to 
the Cauchy formula 

II(i)P(*) dt = - -ImH(t)P(t) dt (2) 

JO 7T 

where the kernel P(t) is an arbitrary analytic function. The integral on the 
left hand side (l.h.s.) is over a circle of radius R. If R is taken large enough, 
we can replace H(t) on the l.h.s. by it's QCD and operator product expansion 
(OPE) counterpart IlQCD(t)- The right hand side (r.h.s.) involves, apart from 
the nucleon pole, an integral over the cut, consisting of a background plus a set 
of nucleonic resonances. Duality means, that the OPE result on the l.h.s. of Eq. 

is equated to the hadronic contribution on the right hand side. Traditionally 
the integrand on the r.h.s. is approximated by the "pole plus continuum" model, 

- Imn(g 2 ) = \ 2 N 6(q 2 - m 2 N ) + -9(q 2 - W 2 ) Imn OPE (g 2 ) (3) 
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Here mjv is the position of the lowest lying pole with residue Xn, the coupling 
of the current to the nucleon state 

(0\r]\n) = A* , 

and an effective continuum threshold W 2 which is determined in the calculation 
and on which the results depend sensitively. 

Most sum rule studies of baryonic currents invoke a Borel transform of the 
correlator, i.e. they use a kernel 

P(t) = e-*/ M2 

which introduces another, more or less arbitrary, parameter providing exponen- 
tial damping of the continuum(when it is small) and suppressing high dimen- 
sional vacuum condensates (when it is large). Stability has to be established 
under variations of the latter parameters. 

The model of Eq.([3]) is rather unrealistic in cases with distinct higher res- 
onances. One may consider the resonances explicitly, e.g. for the nucleonic 
correlator Eq.(UJ one can use in the sharp mass approximation 

- ImriiG? 2 ) = A^r S(q 2 - m %) + y2 A?% 2 - m 2 ) 

7T * — ' 

i 

- Imn 2 (g 2 ) = m N \ 2 N 6(q 2 - m 2 N ) + V mi X^ )2 5(q 2 - m 2 ) - V rm^&tf - m 2 ) 

For the nucleonic correlator considered here, there are four radial recurrences 
iV+(1440), TV" (1535), TV" (1650) and the 7V+(1710), where + or- refers to 
the parity of the state [8] and where the couplings are unknown. The masses 
and widths of these states are well known experimentally.lt is seen from the 
above that the Borel sum rules involving H 2 (q 2 ) are more reliable than the ones 
involving ITi(q 2 ) because of the cancellation between resonances of opposite 
parities |10) . With so many parameters to vary, the sum rule approach is often 
viewed more as an art rather than a science. 

To overcome these intrinsic ambiguities we have introduced some time ago 
a sum rule method [7], originally called ACD, which exploits the analyticity 
properties of the correlator to significantly reduce, in some cases practically 
eliminate, the contribution of the continuum. The breakthrough in the treat- 
ment of the continuum has been the introduction of an integration kernel in 
the FESR tuned to suppress substantially the resonance energy region above 
the ground state. This approach, specially adapted to eliminate pronounced 
resonances, has been recently used to extract very precise values of the light 
quark masses [TT] and condensates [T2]. Our approach is based on the fact, that 
the contribution of the continuum in the integral on the r.h.s. of Eq.([2]) arises 
mostly from the interval 

2.0GeV 2 < t < 3.0Ge^ 2 (4) 
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where the four nucleon resonances lie. This prompts us to choose 

where to = 2.52GeV 2 is the location of the common midpoint of the parity even 
and parity odd resonances. The value of P(t) is small in theregion of the res- 
onances, moreover the contributions of the two parity even resonances iV(1440) 
and jV(1770) (or the two parity odd ones iV(1535) and iV(1650)) come with 
opposite signs, so they tend to cancel in sum rules involving both IIi(g 2 )and 
Tl2(q ). A residual model dependency is still unavoidable as inelasticity, non- 
resonant background and resonance interference are impossible to guess realis- 
tically. Also here our approach helps, as a constant background is eliminated by 
the integration kernel. Having thus minimized the contribution we will neglect 
it. 

The theoretical side of the sum rule, in contrast, is in better shape. The 
correlator (fTJ) is known including radiative corrections and OPE terms up to 
dimension d = 9 [I], [I], [5]. There exists the usual uncertainty about the precise 
values of the condensates and the validity of the factorization assumption used 
for the higher dimensional condensates. We note that the kernel Eq.© will 
introduce only low dimension condensates into the calculation which are known 
(or at least estimated). 

Apart from the references cited above there are a few more attempts to 
evaluate the sum rule ((2J involving a high sophistication on the theoretical 
side which is not always commensurate with the primitive model ansatz on the 
phenomenological side. We therefore think it necessary to present for once an 
(almost) model independent investigation of the nucleonic sum rule. 



2 The calculation 

The invariant amplitudes have poles at t = m 2 



ni(t) 
n a (t) 



X 2 



2 ) 



(t — III 

-m N .X 2 N 

(t - ml) 



Having essentially eliminated the contribution of the continuum, it follows 
then from Cauchy's theorem that 
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|A| 2 P(^) = j U? CD P(t) 

\t\=R 

m N \X\ 2 P{ml) = ^ j n% CD P(t) 

\s\=R 



with 



(27r) 4 nf CD (t) = A t 2 In -4 + A 01 t 2 (in 4 ) + A 4 In -± 

+ ^ 6 - + +A 61 -ln— +A 8 - + .... (6) 

(2^) 4 nf^(t) = Safin ^ + B 7 ~ + 5 9 i + ... (7) 

The coefficients Ai and -Bi are defined as in [5] but for a factor (27r) 4 and 
powers of t taken explicitly 



and 



A 


= — A ^ + -^o),A 01 -- 


A 4 


= -^(aGG) 


A 6 


= (qqqq) 


A 8 


= -^-(27r) 4 Mo (<M«) 


B 3 


= 4tt 2 (qq) (l+|o) 


B 7 


4tt 4 

= - — (ft) (aGG) 


B 9 


= (2, f™a((qqf) 


where a = 


-■ Qs(m2) . The terms B 



-f 



mation. In we have taken 

(Q\qqqaG^a^q\Q) = $ (qqqq) 

with the parameter /i 2 , = 0.8Gey 2 as advocated in [13]. We use the vacuum 
dominance approximation, also for the four-quark condensate, (qqqq) ~ (qq) ■ 
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The question at which scale this relation holds is resolved by allowing for gen- 
erous errors. To avoid the double counting, we keep the logarithmic ln(— <//i 2 ) 
contribution in the polarization operator but neglect its anomalous dimension. 
In any case anamalous dimesion effects are very small [5] . 

For the finite energy sum rule Eq.([2|), we need the known integrals of the 
form 

i ik = ^-j> dttt (H-t)) k 

These are given in convenient form in |14j . 
With our choice of P(t), we then get 



AaR(1-— )- A 6 -A 61 (In _ -_) + -! (8) 



(27r) 4 |A| 2 m„P(r<) = -B 3 R 2 ^- — j B 7 + (9) 

In order to allow the cancellations introduced by P(t) to take place, we choose 
R = 2.92 GeV 2 , i.e. right after the fourth resonance. The nucleon mass is 
obtained by dividing the two sum rules. 



3 Results and concusions 



With standard values of the condensates and using the factorization approxima- 
tion as discussed above, we get a value for the nucleon mass close to the experi- 
mental value. For a more precise statement, we need do discuss possible errors in 
the sum rule. It is important to distinguish two kinds of errors, theoretical and 
experimental. A theoretical error of ±.03GeV arises from the uncertainty in the 
quark condensate. We used (qq) = — (1. 90±0.14) x lO~ 2 GeV 3 at fi = 2GeV cor- 
responding to the limits set by m u + m,d = (8.2±0.6)A/eT^ (at scale fi — 2GeV) 
[TTj in the GMOR relation. Varying the scale parameter /x 2 between AGeV 2 and 
2GeV 2 introduces an additional error of ±0.03GeV. Furthermore there is an 
error due to the strong coupling constant. We use the latest comprehensive 
update analysis at the r-scale [TS] which gives a s (m 2 ) = 0.342 ± 0.01 and leads 
to an error of ±0.6MeV in the nucleon mass. On other theory errors, like the 
one due to the factorization assumption, we have no quantitative handle, but 
under reasonable assumptions they turn out to be less important. We estimate 
the experimental error, by varying to, the zero of the polynomial P(t), between 
2.36GeV 2 and 2.72GeV 2 . The resulting error in the nucleon mass is 

m N = (945 ± 4.5)MeV 

In conclusion, we have presented a simple calculation of the nucleon mass 
using a kernel in the dispersion integral tailored to minimize the contribution 
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of the unknown continuum without involving the higher order unknown con- 
densates. Using standard values of the known condensates, we obtain for m^r a 
value which agrees quite well with the experimental one. 
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